Generation of long-living entanglement between two separate atoms 
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A scheme for non-conditional generation of long-living maximally entangled states between two 
spatially well separated atoms is proposed. In the scheme, A-type atoms pass a resonator-like equip- 
ment of dispersing and absorbing macroscopic bodies giving rise to body-assisted electromagnetic 
field resonances of well-defined heights and widths. Strong atom-field coupling is combined with 
weak atom-field coupling to realize entanglement transfer from the dipole-allowed transitions to 
the dipole-forbidden transitions, thereby the entanglement being preserved when the atoms depart 
from the bodies and from each other. The theory is applied to the case of the atoms passing by a 
microsphere. 
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I. INTRODUCTION 

Generation of entanglement in atomic systems has 
been a subject of intense theoretical and experimental 
study motivated by both the fundamental issue and po- 
tential applications in quantum information processing. 
In this context, the realization of easily controllable long- 
living entangled states of spatially well separated atoms 
has been one of the crucial and challenging problems. A 
number of methods of entanglement preparation between 
atoms have been considered such as the use of quantum- 
correlated light fields interacting with separate atoms, 
thereby transferring their entanglement to the atoms [1- 
6] , appropriate measurements on the light in multi-atom- 
light interaction processes, thereby conditionally project- 
ing the atoms in entangled states [7-17], and the tech- 
nique of quantum reservoir engineering in a cascaded 
cavity-QED setting [18]. 

Photon exchange between two atoms is one of the sim- 
plest processes to entangle two atoms in a common elec- 
tromagnetic field. The effect, which is very weak in free 
space, can be enhanced significantly when the atoms are 
in a cavity [19-21]. Usually attempts are made to mini- 
mize the effect of spontaneous emission. Quite counter- 
intuitively, in certain situations one can take advantage 
of the spontaneous emission for entanglement generation 
[22-25]. Consider, for example, two two-level atoms lo- 
cated in free space with one of them being initially ex- 
cited. This product state is a superposition of a symmet- 
ric (superradiant) state and an antisymmetric (subradi- 
ant) state. If the two atoms are separated by distances 
much smaller than the wavelength, the symmetric state 
decays must faster than the antisymmetric one, leaving 
the system in a mixture of the ground state and the en- 
tangled antisymmetric state. 

The scheme also works at distances much larger 
than the wavelength, if a resonator-like equipment is 
used which sufficiently enhances the atom-field coupling, 
thereby ensuring that a photon emitted in the process 



of resonant photon exchange, which is mediated by real 
photon emission and absorption, is accessible to the two 
atoms. This condition can be satisfied, for instance, when 
the atoms pass by a dielectric microsphere at diametri- 
cally opposite positions [22]. If the distance of the atoms 
from the surface of the sphere becomes sufficiently small, 
then the excitation of surface-guided (SG) and whisper- 
ing gallery (WG) waves can give rise to strong collective 
effects, which are necessarily required to generate sub- 
stantial entanglement. Needless to say that other than 
spherically symmetric bodies can also be used to realize 
a noticeable mutual coupling of the atoms. 

A drawback of the use of two-level-type atoms is that 
the entanglement is transient. In particular, when two 
atoms that have become entangled between each other 
near a body such as a microsphere move away from it 
(and from each other) , then they undergo ordinary spon- 
taneous emission (in free space) , which destroys the quan- 
tum coherence. Preservation of the atomic entanglement 
over long distances between the atoms is therefore not 
possible in this way. 

The contradicting effects of entanglement creation and 
destruction typical of two-level atoms can be combined in 
a more refined scheme involving two three-level atoms of 
A type each (Fig. 1), where the two lower lying states |1) 
and 1 2) such as the ground state and a metastable state 
or two metastable states represent the qubits that are de- 
sired to be entangled with each other [26]. Whereas the 
transition |1) <-> |3) is strongly coupled to the field, the 
transition |2) <-> |3) is only weakly coupled to the field. 
Each atom is initially in the state |1), while the field is 
prepared in a single-photon state. Let us assume that due 
to Rabi oscillations the state |3) of one of the two atoms, 
we do not know which one, is populated. Irreversible de- 
cay to the state |2) is then accompanied with an entan- 
glement transfer forming a (quasi-)stationary entangled 
state between the two atoms with respect to the states 
|1) and 1 2). Its lifetime is limited only by the lifetime 
of the metastable states, and the degree of entanglement 
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achievable can approach 100% in principle. Moreover, 
the scheme is non-conditional and realizable by means of 
current experimental techniques. 

In fact, the model Hamiltonian used in Ref. [26] is 
based on a Dicke-type system and does not allow for 
atoms that are spatially well separated from each other, 
with the interatomic distance being much larger than the 
characteristic wavelengths. However, for many applica- 
tions in quantum information processing or for testing 
Bell's inequalities, large interatomic distances and thus 
the possibility of individual manipulation of the atoms 
are necessary prerequisites. The aim of the present pa- 
per is to close this loophole, by considering two spatially 
well separated A-type three-level atoms appropriately po- 
sitioned with respect to macroscopic bodies, so that the 
two key ingredients enhanced atom-field coupling and 
sharp field resonances can be realized. Note that the sec- 
ond ingredient is absent in the case of a super-lens geom- 
etry [25] . To illustrate the theory, we apply it to the case 
of the two atoms being near a realistic dielectric micro- 
sphere. The formalism used is based on the quantization 
of the macroscopic electromagnetic field and allows to 
take into account material dispersion and absorption in 
a quantum-mechanically consistent manner. 

The paper is organized as follows. In Sec. II the basic 
equations for describing the interaction of N multilevel 
atoms with the electromagnetic field in the presence of 
dispersing and absorbing macroscopic bodies are given. 
In Sec. Ill the theory is applied to the problem of forma- 
tion of an entangled state between two A-type three-level 
atoms. Section IV presents the results obtained for the 
case when the two atoms are at diametrically opposite 
positions outside a microsphere. Finally, a summary and 
some concluding remarks are given in Sec. V. 

II. MASTER EQUATION 

Consider N multilevel atoms at given positions that 
interact with the electromagnetic field in the presence 
of some macroscopic, linear bodies, which are allowed 
to be both dispersing and absorbing. In electric dipolc 
approximation, the overall system can be described by 
the multipolar-coupling Hamiltonian [27] , 

H = Jd 3 rj dwfcjft(r,w)f(r,c«;) 
+ EE &JAm ^ Amm_ E / dw [ d AE( r A,w) + H.c.]. 

Am A 

(1) 

Here, the bosonic fields f(r, u>) and F(r, u), 

[/ fe (r,^),/t,(r>')] = «WJ(w-w , )*(r-r , ) ) (2) 

are the canonically conjugated variables of the system, 
which consists of the electromagnetic field and the bodies 



(including the dissipative system responsible for absorp- 
tion), the RAmn are the atomic (flip) operators 

B-Amn = \m) A A(n\, (3) 

with \m) A being the mth energy eigenstate of the Ath 
atom (of energy hujAm), and 

dA = ^ d Amn R Amn (4) 

m,n 

are the electric dipole operators of the atoms (d^ mn = 
A{m\dA\n) a)- Further, the body- assisted electric field in 
the u> domain, E(r, u>) , expressed in terms of the funda- 
mental variables f (r, u) reads 

E(r,w) = J dVG(r,r>)f(r>), (5) 

where 

G(r, r', w) = iJjL^y/Jmep, u ) G (r, r', w), (6) 

V 7T£0 & 

with G(r, r',w) being the classical Green tensor which 
satisfies the equation 

Vx VxG(r,r»-- e(r,w)G(r,r',w) = «5(r-r') (7) 

together with the boundary conditions at infinity [S(r), 
dyadic S function] . Throughout the paper we restrict our 
attention to dielectric bodies, which are described by a 
spatially varying complex permittivity e(r, u>) =Re e(r, lo) 
+ ilme(r, uj). 

Next we assume that the macroscopic bodies, say, mi- 
crospheres or photonic crystals, act like resonator-like 
equipments such that the excitation spectrum of the 
body-assisted electromagnetic-field shows a resonance 
structure, with the lines being well separated from each 
other. With regard to the atom-field coupling, we as- 
sume that a few atomic transitions can be strongly cou- 
pled to field resonances tuned to them, while all other 
transitions are weakly coupled to the field. Following 
Ref. [28], we decompose the body-assisted electromag- 
netic field into the part (denoted by J '°° dcu . . .) that 
can be strongly coupled to atomic transitions and the 
rest (denoted by J Q du> . . .), which only gives rise to a 
weak atom-field coupling. The Heiscnberg equation of 
motion for an arbitrary operator O that belongs to the 
system consisting of the atoms and the part of the body- 
assisted electromagnetic field that strongly interacts with 
the atoms can then be written in the form of 

pf/OC 

+ lE / AbJ \[d,d A ]MrA,u) 

A L 

+ ^(r A ,uj)[6,d A ]}, (8) 
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where 



//>/oo 
d 3 r J duHw&(T,w)f(T,w) 

+^^fcAm^mra-^ / dui [d A F,(r A , w) + H.C.] . 

A ™ A J 



(9) 



To handle the weak atom-field interaction, i.e., the 
integral J Q °° du; . . . in Eq. (8), we first formally solve the 
Heisenberg equation of motion 

i(r,u;) = - l -[f(v,uj),H} 

= -iwf(r,w) + |^d^G*(r A ,r,o;), (10) 



which yields 

f(r,w,t) = ff ree (r,u;,f) 

+ iE / di'd^OG^r^e-^'-*'), (11) 
ft A Jo 

where ff rco (r, i) evolves freely, 

ffree(r,w,t) = ff rC c(r,w,0)e- Mt . (12) 
Inserting Eq. (11) into Eq. (5), we derive 
E(r,w,t) = E free (r,w,i) 

+ -TEf dt ' e " M)lmG ( r ' r ^ w )dA(*'), (13) 
7re c 2 

where E froo (r, u>, t) is defined according to Eq. (5) with 
ff rco (r, u, t) in place of f (r, u>, t). Introducing slowly vary- 
ing atomic operators 



RAmn(t) — RAmn(t)e lu)Amnt , 



(14) 
(15) 



where the hw An are the atomic energy levels including 
the anticipated media-induced shifts, we may write the 
electric dipole operator, Eq. (4), as 



d-A{t) = d Amn R Amn (t)e 



(16) 



We now insert Eq. (13) together with Eq. (16) in the 
integral J Q °° dcj . . . in Eq. (8), apply the Markov approxi- 
mation to the slowly varying atomic variables in the time 
integral, and take the expectation value. Assuming that 
the free field is initially in the vacuum state, we derive 



(cf. App. A of Ref. [28]) 



(6) = -±([d,H B ]) 



+ i ^2 El (j*AA' ( [A ^Amn] RA'nm) 



A, A' m,n 



+ 

^ AA' {RA'nm \0 , R-An 

~ 2 X] ^2 ( r AA'([d,RA m n]RA>n m ) 
A, A' m,n 

~ ^AA' (RA'nm [O, R-Amn] )) , (17) 

where the primed sum ^2' A A , indicates that A ^ A' and 

the primed sum Yl'm n indicates that transitions that can 
strongly interact with the body-assisted electromagnetic 

field are excluded. In Eq. (17), Hs is defined according 
to Eq. (9), with uj Am being replaced by 



&Am — U Am — S AA , 



where 



rm \ smn 

°AA ~ 2-^i AA> 

n 

with 5 AA being obtained from 

1 f°° 

= - ? V / dujuj 2 

hire c 2 J 



(18) 



(19) 



rmn 
AA' 



d Amn Im G(r A , r A , , uj) d* A , r 

U) — UJ A i mn 



(20) 



(V, principal part) for A=A'. For A ^ A', the parameters 
(5™™, are the dipole-dipole coupling strengths between 
different atoms A and A' . Further, the decay rates T™", 
are defined according to 



_ 2uJ A'mn 

1 heoc 2 



pmn _ 



B(wA'mn) 



x d Amn ImG(r A ,r A ',ui A ' mn )d A , 



(21) 



[O(x), unit step function]. 
Using the relationship 



(O(t))=Tv[p(0)O(t)] 

= Tr[p(t)O(0)]=Tv[g(t)O(0)], (22) 

where p is the density operator of the overall system, and 
g is the (reduced) density operator of the system under 
consideration, and making use of the cyclic properties of 
the trace, from Eq. (17) we derive the following equa- 
tion of motion for the system density operator in the 
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Schrodinger picture: 



nmU 

A, A' m,n 

— RA'nmQRArnn) + H.C. 



^AA' (R-AmnRA'nmQ 



A, A' m.n 

RA'nmQRA run ) 



H.c. 



(23) 



Equation (23) is a generalization of the two-level-atom 
result in Ref. [28] to the case of multilevel atoms. In 
particular, if the conditions 



cmn 
AA' 



°A' A' 



pmn -pmn 

1 AA' ~ 1 A'A 



(24) 
(25) 



are fulfilled, which is the case when, for example, the 
atoms are identical and located in free space or at equiv- 
alent positions with respect to the macroscopic bodies, 
then the master equation (23) takes the somewhat sim- 
pler form of 



Q = - 



^ + Hr> , g 

2 ^AA'* {RAmnRA'nmQ 

A, A' m.n 

2RA'nrnQRAmn + QRAmnRA'nm), (26) 



where 



#D — — ^2 X/ ^AA' RAmnR-A'r 



(27) 



A, A' m>n 



describes the dipole-dipolc interaction between the 
atoms, with A™^, being the dipole-dipole coupling 
strengths, 



A mn 
^AA' 



cmn I x nr ' 
AA' + A' 



A' A- 



(28) 



According to Eq. (26), the (undamped) system is gov- 
erned by an effective Hamiltonian equal to H§ + Hd. 
Note that this is not true in general, but only under the 
conditions (24) and (25). 

To construct the (formal) solution to the master equa- 
tion (26), we first rewrite it in the form of 



g = Lg + Sg, 



(29) 



where L and S are superoperators which act on g accord- 
ing to the rules 



Lq=--(Hq-$W), 



Sg — ^AA' RA'nmQRArnn, 



(30) 
(31) 



A, A' m,n 



and the non-Hermitian "Hamiltonian" TL reads 

H =H§+ Hb - y ^2 X/ ^AA' RAmnR-A'nm- (32) 
A, A' m>n 

From Eqs. (29)-(31) it then follows that 

g(t) = e^-^gito) + [ dtie 1 ^*^ S g{ti). (33) 
By iteration, from Eq. (33) one readily finds 

oo 

£(t) = X> (n) (*), ( 34 ) 



n=0 



where 

0(°>(t) = e £ (*-«°>0(t o ), 



(35) 



Jto Jto J to 



x 5'e i( *"- t "- l) . . . 5e i(tl -* o) ^(t ), n = 1, 2, 3 . . . . 

(36) 

Although Eq. (34) is not a perturbative expansion, it can 
be helpful, as we shall see below, in finding the explicit 
solutions to the master equation. 



III. TWO THREE-LEVEL ATOMS OF A TYPE 

A. Solution to the master equation 

Let us specify the atomic system and consider two 
identical three-level atoms A and B of A type as sketched 
in Fig. 1. We assume that the dipole-allowed tran- 



13} 




12} 



11} 



FIG. 1: Three-level atom of A type. 

sition |1) <-» |3) is tuned to a well pronounced body- 
induced electromagnetic field resonance, thereby giving 
rise to a strong dipole-allowed atom-field coupling. Fur- 
ther, the dipole-allowed transition |2) <-> |3) is assumed to 
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be weakly coupled to the body-assisted electromagnetic 
field, and the transition between the states |1) and |2) is 
dipole-forbidden. Restricting our attention to two atoms 
at equivalent positions with respect to the macroscopic 
bodies, so that corresponding transition frequencies are 
equally shifted and the relations 



a 31 _ A 31 a 32 _ A 32 



(37) 



p32 _ T^32 p32 _ t^32 /oo\ 

1 AA 1 SB' 1 AB — 1 BA l d0 J 

hold [cf. Eqs. (24) and (25)], we may apply the master 
equation in the form of Eq. (26) and its solution in the 
form of Eqs. (34)-(36), with Eqs. (31) and (32) being 
explicitly given by 



Sg= ^2 F A 2 , A „R A „ 23 qRa> 32 



(39) 



and 



A'A"=A,B 



lk 



H =Hs + H D — F a 2 , a „Ra>32Ra"23, (40) 



A' A' 1 



where 

H s = J d 3 r J dwfiwf+(r,w)f(r,o;) 

3 

mm 



A'=A,B m=l 
/•foe 



p/OO 

- ^2 dw[dA'3i-RA'3iE(rA',w) + H.c], 

A'=A,B ® 



Hd = —(hA AB R A 3iRB13 + A j4B i? J 432i?B23) + H 



(41) 



(42) 



with the rotating-wave approximation in Eq. (41). 

To specify the initial condition at time to, let us assume 
that the two atoms are initially in the ground state |1, 1) 



(\i,j) = \i) A <E>\j) B , i,j = l,2,3) and the rest of the system 
is prepared in a state 



\F) 



= J duj J d 3 rC(r,w,t )f t (r,w)|{0}>, (43) 



where C(r, u>, to) as a function of ui is non-zero in a small 
interval around u> ~ u)A3i = &B31, and |{0}) is vacuum 
state with respect to this frequency interval. The initial 
density operator can then be given in the form of (to = 0) 

0(0) = |*(0))<*(0)|, |*(0)) = |1,1)®|F). (44) 

In order to determine the density operator at time t, we 
begin by calculating the first term of the series (34), viz. 

g^(t) = e Lt g(0) = Mt))m)l (45) 
where the (damped) state vector 



obviously obeys the equation 



(46) 



(47) 



Recalling the initial condition (44) and the form of Tt, 
Eqs. (40)-(42), we may expand |*(t)) as 

|*(t)) = CsiWe-^+^tfi, 1) ® |{0}) 
+ C 13 (t)e- i ( QA3+QB1 »\l,3)® |{0}) 
+ J duo J ^ re -i(Q A1 +u, B1 +u,)t 

x C(r,w,i)?t(r,w)|{0})<8>|l,l). (48) 

We now substitute Eq. (48) into Eq. (47) and make ex- 
plicitly use of Eqs. (40)-(42). Straightforward calcula- 
tion yields the following system of differential equations 
for the expansion coefficients: 



C31 = -h T AAC 3 i +iA A 1 B C 13 + l -J^ dejj d 3 rd A 3iG(r A ,r, W )C(r, W ) e -^-^)*, (49) 
C w = -\T^ B C w + iti% B C zl + ±J dcu yd 3 rd B3 iG(r B ,r,^)C(r,^) e -^-^ 31 )*, (50) 

C(r, w) = l e *(«-^si)* [d ^ 31 G*(r A , r, W )C 31 + d^ 31 G*(r B , r, W )C 13 ]. (51) 

I 

Recall that &A31 =&B3i- Inserting the formal solution differential equations 
to Eq. (51) in Eqs. (49) and (50), we derive, on making 

use of the properties of the Green tensor, the integro- C31 = — 2^ aa^S! + iA AB C\3 

+ / dt'[K AA (t-t')C3i{t') 
Jo 



+ K AB (t-t')C 13 (t')]+F 31 (t), (52) 



6 



C13 - - ^r 3 z B Ci 3 + iA 3 AB Q 



AB°31 

+ / dt'[K BB (t-t')C 13 (t') 
Jo 

+K BA (t-t')C sl {t')]+F ls {t), (53) 
where the kernel function Ka'A" (t) is defined by 

K A , A „{t) = -- / do; — e -*(«-^3i)* 

x cU'3i ImG(r j4 ',r A //, w) d^„ 31 (54) 

[A'(A") = A, B], and the free-field driving terms F 31 and 
F13 read 

F 3 i(t) = % -J^ du;j d 3 rd A31 G(r A ,r, W ) 

x C(r,w,0)e- l(w -^ 3l) *, (55) 

Fia(t) = \j Q duj J d 3 rd B3 iG(rs,r,w) 

x C(r,oj,0)e- l(ul - QB31 '> t . (56) 

Note that for identical atoms at equivalent positions with 
respect to the macroscopic bodies 

K AA (t) = K BB (t), K AB {t) = K BA {t). (57) 

Instead of considering the probability amplitudes C31 and 
C13, it is advantageous to introduce the probability am- 
plitudes 



^ = 2^ (C 31 ± C 13 ) 



(58) 



which are the expansion coefficients of \^>) with respect 
to the atomic basis 



|±i 3 ) = 2-s(|3,l>±|l,3», 
so that Eq. (48) takes the form of 

|*(t)) = Cf^e-^+^'l+is) ® |{0}) 
+ C^i^e-^+^^l-u) ® |{0}) 



(59) 



j3 r e -i{£o Bl+^)t 



xC(r,w,t)ft(r,w)|{0})®|l,l>. (60) 

From Eqs. (52)-(58) it is not difficult to see that the 
differential equations for C_J_ 3 decouple 

A13 _ / _i_ - A 31 l r 32 \ r 13 

+ f dt'K±(t-t')Cl 3 (t') + F±(t), (61) 
Jo 



where 

Jf±(t) = JfAA(t)±«- AB (t), (62) 
F±(i)-2- 1 / 2 [F 31 (t)±F 13 (t)]. (63) 

The field resonance strongly coupled to the atomic 
transition |1)<-»|3) can be typically modeled by a 
Lorentzian, with ujq ~ u>A3i and Aojq being the central 
frequency and the half width at half maximum, respec- 
tively. In this case, Eq. (54) can be approximated by 

K A , A „(t)=-T 3 A 1 , A „e- i ^- QA ^ t 



2tt/ 



duj 7 ^3 t t - . 64 

(uj - w c 2 + Aw£ 



where T 3 A , AII is defined according to Eq. (21), but with 
u) A si being replaced by u>c, 

2uj 2 

T^A'A" = j-^d A , 3 ilmG(r A ,,r A ,,,cjc)d A ,, 31 . (65) 
From Eq. (64) it then follows that (t > 0) 

K A >A»(t) = -|r 3 4 1 , A „Ao;ce- < ( A - <Awc ) t (66) 

(A = u>c — & A 3i). Using Eq. (66) and differentiating both 
sides of Eq. (61) with respect to time, we find that C]_ 3 
satisfies the second-order differential equation 

C 13 + a 1± C 13 + a 2 ±C 13 = F±{t) + i(A - iAoj c )F±{t), 

(67) 



where 



oi± = *(A T A AB ) + Alo c + ir 3 ^, 



(68) 



a 2 ± = g± + (A - iAu; c )(±A AB + i\T 3 A \), (69) 



with 



2 1t^31a,, t^31 r3i 4- r3i 
.9± = 2 1 ± Acj C, 1± =l AA ±l AB - 



(70) 



If C±(t) are known, then the probability amplitude 
C(r, u>,t) can be obtained from Eq. (51) together with 
Eq. (58). 

To calculate the terms g( n \t) (n > 0), Eq. (36), of 
the series (34), we note that the action of the opera- 
tor S, Eq. (39), on = corresponds 
to atomic transitions |3) — > |2). Thus, only the states 
1 1 , 3) and |3, 1), or equivalently |±i 3 ), can contribute to 
S[\^(t))(^(t)\]. It is not difficult to see that 

S(|± 13 }(± 13 |) = 

r 3 4 2 4 |± 12 )(± 12 | T |r 3 _ 2 (|+ 12 )(+ 12 | - 1-12X-12I), (7i) 

5(|± 13 )(Ti 3 |) = 

T 32 A \± 12 )( Tl2 \ - |r 32 (|± 12 )( T i2| - |Ti2><±i 2 |) (72) 
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[rf = T% A ±Tf B , |±i 2 ) = 2-i(|2,l)±|l,2))]. Combining 
Eqs. (60), (71), and (72), we derive 

Sg^(t) = S[\^(t))m)\] 

= |{0})({0}|®{(irf|c| 3 | 2 

+ l r 3_ 2|c l 3|2)| + i2)(+i2| + ( l r3+ 2 |c 13 |2 

+ iri 2 |c| 3 | 2 ) I-12X-12I + [^rfcfc! 3 * 

+±r 32 C 1 + 3 *C 13 ) I+12X-12I +H.c.]}, (73) 

SS(\H>(t))(*(t)\)=0. (74) 

Recalling that H, Eqs. (40)-(42), acts on atomic states 
in the subspace spanned by |±i 3 ), we see that 

e L ^S(\*(t 1 ))(*(t 1 )\)=S(\*(t 1 ))(*(t 1 )\), (75) 
leading to 

e W(t)= [* dtjimwmtji) (76) 

Jo 

[cf. Eq. (36)]. Further, Eqs. (74) and (75) imply that 
g( n ) =o if n > 2. Thus, the solution to the master equa- 
tion reads 

m = i*(t)x*(t)i + f *tiS\mi))mi)\] (77) 

Jo 

together with Eqs. (60) and (73). 

B. Stationary limit 

Let us restrict our attention to the stationary limit 
t— > 00. Since i 7 3i(t) and -Fi3(t) approach zero as t tends 
to infinity, Eqs. (49) and (50) imply that 

lim C^(t) = 0. (78) 

t— ^00 

Inserting Eq. (60) in Eq. (77) and taking the trace with 
respect to the f -field, we derive 

Tr fie i d £(t-> 00) = g at , (79) 

Qat = «+|+i2) (+12I + a-|-i2>(-i2| 

+ (/JI+12X-12I + H - c -) + (!"«+ - 1|, (80) 

where 

roo 

a ± = dt(ir^|Ci 3 | 2 + ^|Ci 3 | 2 ), (81) 
Jo 

POO 

0= dt (irfcf cl 3 * + §ri 2 c* 3 *<7l 3 ) . (82) 

Jo 



To determine the accessible entanglement of the two 
atoms, it may be instructive to study the concurrence of 
the atomic subsystem, which may be regarded as being a 
measure of entanglement [29]. For this purpose, we have 
to calculate the spin-flipped density operator 

Pat = {0Ay ® a By ) Pit (°By ® ^ Ay) , (83) 

where 

(°A { B)y) mn = ( i "J ) ( 84 ) 

[m(n) = 1,2], and to determine the two nonzero eigenval- 
ues A± of PatPat- A somewhat lengthy but straightfor- 
ward calculation yields 

A± = \ {a 2 + + a 2 _-2 [(Re/?) 2 - (Im/3) 2 ] } 

± Wi( a + + a -) 2 - 4(Re/3) 2 ] [(a+-a_) 2 + 4(Im/3) 2 ] , 

(85) 

which then determine the concurrence 

C = ^/\^-^/\Z, (86) 

the value of which is in the interval [0,1]. The nearer 
1 the value of C is, the higher is the degree of entangle- 
ment. Equations (85) and (86) reveal that a noticeably 
entangled state of the two atoms can be generated if 

a+(a_)> a-(a+),\(3\, (87) 

thus C — ► a+(a_). Needless to say that the entangle- 
ment condition (87) is already expected from inspection 
of Eq. (80). 

C. Different Coupling Regimes 

Let us return to Eq. (67) and focus on the case where 
F± (t) ~ -i(A - iAcj c )F± (t) (88) 

is valid, so that the term on the right-hand side in 
Eq. (67) can be omitted. Obviously, this is the case when 
initially the (Lorentzian) field resonance of mid-frequency 
u>c and width Aloc is excited (for details, see Sec. HID). 
Under the initial conditions 

Ci 3 (0) = 0, C 1 ± 3 (0)=F ± (0), (89) 

the solution to Eq. (67) can then be written in the form 
of 

C ™(t) = ^±iPi e -« 1± */2( e9± t/2 _ e -9±t/2) ) ( 90 ) 

Q± 

where 



1± = Ja 2 1± -4a 2 ± • (91) 
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Restricting again our attention to the stationary limit, 
we further assume, for simplicity, both the detuning A 
and the dipole-dipolc coupling strength A Z AB vanish, i.e., 
A = and A^g = 0. Since even under these conditions the 
explicit form of the expansion coefficients a±, Eq. (81), 
and f3, Eq. (82), is rather involved, we renounce its pre- 
sentation here but consider instead some instructive spe- 
cial cases. 

From Eqs. (68) and (90) it is seen that the damping 
constant of C_|_ 3 is determined by the sum of the half 
width at half maximum of the field resonance strongly 
coupled to the transition |3) «-> |1) and the half width 
at half maximum of the transition |3) — > |2), Aloq and 
T A 2 A /2, respectively. Due to the finite Ao>c, an atom 
tends to occupy the state |1), while the effect of the finite 
T A 2 A is that the atom prefers to occupy the state |2). We 
may therefore restrict ourselves to situations in which 



(92) 



To achieve noticeable entanglement, the interatomic cou- 
pling should be sufficiently strong, i.e., \T A X B 
\ t ab I -^TaA) equivalents, 



T 3 aa and 



r^ 1 

1 T 



> 1, 



L ±(T) 
p32 



> 1. 



(93) 



Note that the first inequality is equivalent to g± » g T 
[cf. Eq. (70)]. We now distinguish between the following 
three cases. 



(b) g± > g T > Tf A > Auj c 

When both g± and g^ dominate the other parameters, 
then the states I+13) and | — 13) are both strongly cou- 
pled to the medium-assisted electromagnetic field, and 
Eq. (90) approximates to 



M _ F ±R e -rf A t/4 sin{g±t ^ (9g) 
9± 



which is exactly analogous to Eq. (94) . The steady-state 
density operator parameters a± and (3 take the approxi- 
mate form of 

a+ ir 3 +2 l^ + (0)l 2 + i r32 l^-WI 2 m 

"± — 2 L ± ^2^32 ^ 2 L T r ,2 p32 V^ 3 > 

9+ L AA .9- 1 AA 

and, for # +( _) »#_(+), 



13 = [TfF+ (0)F*(0) +T 32 F* (0)F_(0)] 



r 32 



AA 



(100) 



(-) 



( c ) r 3 A 2 4 » 5± » 5T , Alu c 

When the value of T A 2 A sufficiently exceeds the values 
of the other parameters, then from Eq. (90) it follows 
that 



cm 



2F±(0) 

p32 
1 AA 



e -\u c t _ e -T 3 / A t/2 



(101) 



(a) g± > T\ 2 A > Alo c > g T 

In this case, either the symmetric state I+13) or the 
antisymmetric state |— 13) is strongly coupled to the 
medium-assisted electromagnetic field whereas the other 
one is weakly coupled. For the strongly and weakly- 
coupled states, respectively, Eq. (90) approximates to 



cm 



^e- r -*/ 4 sin( 5± i), 



9± 



and 



cm 



2jy(Q) 

r32 

1 AA 



e -A^ c t _ e -T A \t/2 



(94) 



(95) 



It is seen that C± (t) undergoes damped Rabi oscillations 
of frequency g±, while C^(t) undergoes a two-channel 
exponential decay. The steady-state density operator pa- 
rameters a±, Eq. (81), and /3, Eq. (82), approximate to 



a± 



1 F 32 
2 L ±(T) 



\F H -)(0)\' 
n 2 p32 

#+(-) aa 



|F-( +) (0)| 2 



/3 = [r3_ 2 F + (o)F*(o) + ri 2 iq:(o)F_(o)] 

for >£_(+). 



p32 
1 AA 



2g 4 



(96) 



(97) 



+(-) 



i.e., the behavior typical of weakly-coupled states is ob- 
served [cf. Eq. (95)]. In this approximation, the steady- 
state density operator parameters a± and (3 read 



a± 



-^32 



l^+(0)| : 



± (Tf A r(Au Jc + 2gl/rf A ) 
r32 l^(0)| 2 

+ * (r 3 4 2 j 2 (A^ + 2. 9 2 /r3 2 4 ) 

and, for # +( _) » #_(+), 

TfF + (0)F*(0) + r 32 F*(0)F_(0) 



n 2 4 ) 2 (A^c + .9 2 ( _ ) /r 3 4 2 4 ) 



(102) 



(103) 



D. Preparation of the initial state 

One possible way to initially prepare the medium- 
assisted electromagnetic field in the desired quantum 
state (43), is to use an additional atom, say atom D, 
such that o>_D3i = <l>A3i = &B3i = ^c- Let the transition 
|1) <-» 1 3) of atom D strongly interact with the medium- 
assisted electromagnetic field in the absence of atoms A 
and B. This can be achieved, for instance, by using 
atomic beams and letting atom D pass the equipment 
before atoms A and B pass it. When atom D initially 
prepared in the excited state |3) strongly interacts with 
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the medium-assisted electromagnetic field initially pre- 
pared in the vacuum state, then an interaction time can 
be chosen after which the atomic excitation is transferred 
to the field. 

The probability amplitude of finding, after some inter- 
action time At, atom D (regarded as an effective two- 
level system) in the ground state and the f-field in a 
single-quantum state is [28] 

C(r,cu,t = 0) 

= l - [ At' d* mi G*{v D , r ,u>)S u - a ™V 'C UD (t'), (104) 

n J-At 

where 

C Ud (t) = e -A"c(*+At)/2 cos[gD{t + At)] (105) 

is the probability amplitude of finding the atom in the 
upper state. Here, 



9d = v /r^Ao;c/2 



(106) 



is the single-atom Rabi frequency, with T^ D being deter- 
mined according to Eq. (65). Substitution of Eq. (104) 
into Eqs. (55) and (56) yields 

,o 

F 31 (t)= dt'K AD (t-t')C UD (t'), (107) 

J -At 
,0 

F 13 (t)= dt'K BD (t-t')C UD (t'), (108) 

J -At 

where Kbd^) is defined according to Eq. (54). Note that 
F±(t), Eq. (63), calculated by using F31 and F13 given in 
Eqs. (107) and (108) fulfills Eq. (88). To calculate F±(0), 
we fix the interaction time At such that Cjj d (0) = 0, thus 



At = 



(109) 



Combining Eq. (63) with Eqs. (105)-(109), we derive, 
on applying the Lorentz approximation according to 
Eq. (64), 



F ± (0) = -i=^±exp(-Ac " 
V2 9D V 2 9D 



where 



9D± = \(T^ D ±T^ D )Au:c/2 



(110) 



(111) 



and T^2 D and are defined according to Eq. (65). 

In Eq. (110), the exponential factor characterizes the 
photon loss during the interaction time due to the finite 
width of the field resonance. Obviously, the better the 
strong-coupling condition Awc<Cgu is fulfilled, the less 
the photon loss is. In particular, when atom A (or B) 
changes places with atom D and the orientations of the 



transition dipole moments of atoms A (or B) and D are 
the same, then from Eq. (110) it follows that (Au>c <C go) 



F±(0)*-g±, F T (0)~-4/ 3± . 



(112) 



It is worth noting that, as we will see in Sec. IV, the 
highest degree of entanglement can be achieved in case 
of equal positions of atoms D and A (or B). 



IV. ATOMIC ENTANGLEMENT NEAR A 
DIELECTRIC MICROSPHERE 

Let us apply the theory to two atoms near a dispers- 
ing and absorbing dielectric microsphere (of radius R) 
characterized by a Drude-Lorentz type permittivity 



e(w) = 1 



(113) 



(ujp, coupling constant; lot, transverse resonance fre- 
quency; 7, absorption parameter), which features a band 



gap in the region u>t < w < lul 
Ree(co) <0. 



\/^\ + u)'p, where 



A. Two-atom coupling 

Making use of the Green tensor for a dielectric sphere 
[30], one can show, on assuming radial dipole orienta- 
tions, that Eq. (21) leads to 



^A'A" — ^A'A" — frpRcy^ 



1=1 



«' + ,yg +i V '(»r) 

(kr) 2 ' 



31 



(kr) + Bf T (u>)h\ 1 \kr)] P;(cos6») (114) 



[u> = u>%r = > 0; k = uj/c; r = r A ,= r A „ (> R), ra- 
dial position of the atoms]. Here, Tq is the single-atom 

decay rate in free space, ji(z) and h\ [ (2) are the spher- 
ical Bessel and Hankcl functions, respectively, Pi(x) is 
the Legendre function, 9 is the angle between the two 
transition dipole moments (Id^'mnl = |d J 4"mn|), and the 
scattering coefficients B^ (w) read [30] 



r(w)ii(«2)[«iji(«i)]' - Ji(zi)[z2ji(z2)}' 



eMjiMlzihPiz!)]' - hf'Wfojtizi)]' 



( 115 ) 

where z-i — kiR, k\ = k, and k^ — y/s(co)oj/c. Note that ra- 
dially oriented dipoles couple only to TM waves, whereas 
tangentially oriented dipoles couple to both TM and TE 
waves (for details, see, e.g., [31]). Needless to say that 
9 = in case of a single atom (A' = A"). 

The complex roots of the denominator of the reflec- 
tion coefficients Bf (uu) determine the positions and the 
widths of the sphere-assisted electromagnetic field reso- 
nances. When lj coincides with a resonance frequency, 
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FIG. 2: The two-atom collective decay rate T^, [Eq. (114), 
A' 7^ A"] as a function of the angle 9 between the transition 
dipole moments for u> = 1.0501 lot- The two atoms are at 
distances Ar = r — _R = 0.14At (\t = 2ttc/ujt) from the surface 
of a dielectric sphere (cup —0.5u>t, 7 = 10~ 6 ojt, R=W Xt)- 



say wc, then the corresponding I term in Eq. (114) is the 
leading one, thus 



1 A' A" — ~ l ' 



oRej 



Z(/ + l)(2/ + l) 



(kr) 



K'{kr) 



ji(kr) + B? \u)h\ 1] '(At) Pz(cos#)S> (116) 



(w ~ wc)- Equation (116) implies that when the two 
atoms (A' ^ A") are at diametrically opposite positions 
with respect to the sphere, i.e., 8 — tt and hence P/(cos#) 
= (— 1)', then the interaction of the symmetric (antisym- 
metric) state with the sphere-assisted electromagnetic 
field is enhanced, while the antisymmetric (symmetric) 
state almost decouples [cf. Eq. (70)]. 

The dependence on 9 of Ta'A" {A' ^ A") as given by 
Eq. (114) is illustrated in Fig. 2, where the atomic tran- 
sition frequency ui is chosen to be close to a microsphere 
resonance frequency. From Figs. 3 and 4 it is clearly 
seen that the value of T + (r_) can drastically exceed 
the value of T_ (r + ) when the two atoms approach the 
microsphere and the transition frequency equals a reso- 
nance frequency. Recall that T + (r_) is a measure of the 
strength of coupling of the symmetric (antisymmetric) 
state to the sphere-assisted field. In particular, Fig. 4 re- 
veals that there is an optimum distance - the distance at 
which the solid curve attains the minimum - for which 
the best contrast between F + and T_ can be realized. 
With increasing distance of the atoms from the sphere, 
the values of both T + and T_ tend to the free-space value 
r as they should. 

Figures 2-4 refer to atomic transition frequencies 
within the band gap. In this case, the strong two-atom 
interaction observed when the atoms are at diametrically 
opposite positions with respect to the sphere is mediated 
by SG waves. Of course, the effect of enhanced T + (r_) 
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FIG. 3: The two-atom decay rates F+ = T A i A i + Ta'A" (solid 
curve) and r_ = T a i a i — Y A 'a" (dotted curve) for the sym- 
metric and antisymmetric states, respectively, as functions of 
the transition frequency w, with Fa'a" from Eq. (114) for 
9 = 7v. The other parameters are the same as in Fig. 2]. 




FIG. 4: The two-atom decay rates F + = T A i A i + T A i A n (solid 
curve) and T_ = T A * A i — T A i A h (dotted curve) for the sym- 
metric and antisymmetric states, respectively, as functions of 
the distance Ar between the atom and the surface of the di- 
electric sphere, with T A i A n from Eq. (114) for 8 — tt. The 
other parameters are the same as in Fig. 2] . 



and simultaneously reduced T_ (r + ) can also be observed 
for transition frequencies below the band gap. In this 
case, the cavity- assisted field resonances correspond to 
WG waves. An example is shown in Fig. 5. Figures 3 
and 5 also convey a feeling of the sharpness of the field 
resonances, which ranges from being very sharp to being 
less so. The sharpness can be improved by increasing 
the microsphere radius or by reducing the material ab- 
sorption. Note that WG waves much more suffer from 
absorption than do SG waves (see, e.g., Ref. [31]). 
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5: The two-atom decay rates r + = T A i A i + T A i A h (solid 
curve) and F_ = T A i A i — T a >a" (dotted curve) for the sym- 
metric and antisymmetric states, respectively, as functions of 
the transition frequency to, with T A 'A" from Eq. (114) for 
=7r. The other parameters are the same as in Fig. 2]. 



B. Entanglement of two A- type atoms 

The results given in Sec. IV A show that the optimal 
positions of two A- type atoms A and B, which are de- 
sired to entangle with each other near a microsphere, are 
diametrically opposite with respect to the sphere. Fur- 
ther, the transition frequency Qa31 = &B3i should coin- 
cide with the (mid-)frequency ujc of a sufficiently sharply 
peaked sphere-assisted field resonance, so that the strong- 
coupling regime is realized and the first of the condi- 
tions (93) is satisfied. Finally, the transition frequency 
&A32 =^B32 should coincide with the (mid-)frequency of 
some moderately peaked sphere-assisted field resonance, 
so that the second of the conditions (93) is also satisfied, 
but the weak-coupling regime applies, thereby giving rise 
to an irreversible decay channel. As a result, the condi- 
tion (92) can also be expected to be satisfied. By choos- 
ing atoms with appropriate transition dipole matrix ele- 
ments well matching cavity-assisted field resonances (for 
more detailed estimations, see Ref. [28]), all the condi- 
tions including both the inequalities characterizing the 
three cases (a)-(c) in Sec. Ill C and the field-preparation 
conditions (112) can be fulfilled. Let us examine the cases 
(a)-(c) in more detail. 



(a) g± > > Auj c > g T 

For defmiteness, let T^ 1 > V 3 } and T^ 2 >T 32 . When 
atom A (or B) changes places with atom D, which pro- 
vides the initial field excitation, and Eq. (112) applies, 
then Eqs. (96) and (97) lead to 



a+ ~ 1, 



p32 



«1, 



(117) 



(118) 



< 1. 



(119) 



Hence, an almost perfectly entangled state is produced, 
£at^ I+12X+12I [sec Eq. (80)], and, accordingly, C~l is 
achieved. Clearly, a+ = 1 (C = 1) cannot be exactly re- 
alized, because of the losses unavoidably associated with 
the always finite width of the field resonance. It is worth 
mentioning that when the positions of atoms D and A 
(or B) are different from each other (e.g., atom D was 
equidistant from atoms A and B), then the degree of 
entanglement that can be achieved is smaller than that 
in case of equal positions in general. Note that when 
r 31 »r 31 and T 32 > rf, then g at ~ |+ 12 )(+ 12 | is also 
valid. For T^ 1 >r 3 F 1 and > T^ 2 , however, the roles 
of a + and a_ are interchanged and g at ~ | — 12) (— 12 1- 

In the scheme, the two-atom system undergoes, e.g., 
fast |1, 1) <-» I+13) Rabi oscillations as long as one of the 
two atoms jumps to state |2), but we do not know which 
one. Hence, the result is the entangled state between 
one atom in the state |2) and the other in the state |1). 
The time after which the stationary limit is established is 
determined by the lifetime ~ (r^) -1 of the short-living 
state I+13), while the long-living state | — 13) of lifetime 
<~ (Aivc)^ 1 is practically unpopulated [cf. Eqs. (94) and 
(95)]. 



(b) g± > g T > Tf A > Aloq 

For defmiteness, we again assume that r 3 ^ 1 ^> r 31 and 
r 32 >r 32 From Eqs ^ ( 99 ) and ( 10Q ) togcther witn 

Eq. (112) we obtain 



r 32 2g 2 

+ 4p < 1, 
9+ 



< 1. 



(120) 



(121) 



(122) 



Thus, this coupling regime leaves the two atoms in an 
entangled state analogous to case (a). However, since the 
inequality g± 3> g^ requires g± to be as large as possible 
and g T to be as small as possible, while the inequality 
g T » T 3 ^ requires that g T is much lager than T A 2 A , it 
may be more difficult to realize this regime. Note that 
g T is the smallest or one of the smallest parameters in 
the cases (a) and (c). 



(c) Tf A >.g± >.g T ,Aw c 

In this case, the irreversible decay from state |3) to 
state 1 2) is so dominant that Rabi oscillations are fully 
suppressed in the time evolution of both C+ 3 and Ci 3 [see 
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Eq. (101)]. From Eq. (102) we obtain, on again assuming 

p31 > p31 and p32 ^ p32 and making use f Eq ( 112 ) ; 



a+ 



2 9 l/rf A 



(123) 



To generate the entangled state I+12), i.e., a + ~ 1, the 
additional condition 



9+ > A ^C 



r 32 



AA 9+ 



(124) 



must be required to be satisfied, as can be seen from 
Eq. (123). The parameters a_ and (3 then read 



p32 2 
' • ^ T^T- + — « 1, 



2r 32 



.9+ 



/3~ -4p < 1. 
9+ 



(125) 



(126) 



In a similar fashion, it can be shown that in case of 
T^ 1 > T^ 1 and T 3 ^ > T^ 2 the antisymmetric entangled 
state |— 12) is generated. 

The inequality (124) can be understood as follows. For 
F+(0)~ -g+, Eq. (101) yields 



C?(t) ~ -(2 5+ /r&) e 



,-rf A t/2 



(127) 



i.e, Ci 3 (i) ~ g+/Y A 2 A . Thus, though one can allow for 
g + /Y AA <§C 1, this ratio has still to satisfy the inequality 
(124) such that there is a nonvanishing probability that 
one of the atoms can reach the state |3) from the initial 
state |1) to jump to the state |2). 



V. SUMMARY AND CONCLUSIONS 

We have proposed a scheme for non-conditional prepa- 
ration of two spatially well separated identical atoms in 
long- living highly entangled states. The scheme uses A- 
type atoms passing a resonator-like equipment of realis- 
tic, dispersing and absorbing macroscopic bodies which 
form electromagnetic field resonances, the heights and 
widths of which are determined by the radiative and non- 
radiative (absorption) losses. The lowest lying atomic 
state and the lower lying excited state, which can be the 
ground state and a metastable state or two metastable 



states, play the role of the basis states of an atomic 
qubit. The atoms initially prepared in the lowest ly- 
ing states, are pumped by a single-excitation "pulse" of 
the body-assisted electromagnetic field, thereby strongly 
driving the dipole-allowed transition between the low- 
est and highest lying atomic states. In this way, one 
of the two atoms - we do not know which one - can 
absorb the single-photonic excitation, and subsequent ir- 
reversible spontaneous decay of the excited atomic state 
to the lower lying excited state, the transition of which 
to the lowest lying state is dipole-forbiddcn, results in a 
metastable two-atom entangled state. 

To be quite general, we have first developed the theory, 
without specifying the atoms and the equipment whose 
body-assisted electromagnetic field is used for the the col- 
lective atom-field interaction. For the case of two A-type 
atoms, we have derived the general solution of the cou- 
pled field-atom evolution equations and presented special 
coupling conditions under which high-degree entangle- 
ment can be achieved. We have then applied the the- 
ory to the problem of entanglement of two A-type atoms 
near a microsphere. In particular, we have shown that 
the scheme is capable of realizing strong coupling in one 
arm and weak coupling in the other arm of the A con- 
figuration. In this context, we have also analyzed the 
preparation of the initial single-photonic field excitation 
required for initiating the process of entanglement. 

In contrary to the common sense that the existence 
of dissipation spoils the quantum coherence of a sys- 
tem, dissipation is here essential to transfer the entangle- 
ment from the strongly driven transitions to the dipolc- 
forbidden transitions. The fact that only ground or 
metastable states serve as basis states of the qubits guar- 
antees the long lifetime of the entangled state. It is worth 
noting that the scheme renders it possible to test nonlo- 
cality for a two-atom system. An atomic pair passing by 
a microsphere and being entangled there, can be sepa- 
rated from each other and one can be sure that in the 
meantime the entanglement is not lost. 
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